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Abstract- In the present investigation a steady-state 
mathematical model is developed for a typical four stage 
CEHH bleaching sequence in Indian Pulp and Paper 
Industry. Each stage of the sequence is composed of 
chemical additions and mixing, reaction in a retention 
tower, and washing. Unit operation models were 
formulated using mass balances on liquor, fibres, kappa 
number, chemicals and COD based on the assumption of 
perfect mixing and quasi steady-state, The COD is used 
in this paper as an indication of the amount of effluent 
coming out from four bleaching stages, COD 
measurement of bleach plant effluent separately is not a 
routine mill measurement; so COD is calculated in the 
basis of correlation of COD given by Anjana (9). In order 
to yalidate the model, the simulated results were 
compared the literature data; a good agreement between 
the simulation results with the literature data is achieved. 
In the case studies, the optimization of the process 
variables are performed in which the effects of operating 
temperature, residence time, consistency, bleaching 
liquor flow rate and the wash water flow rates in the 
washers are studied. 


Keywords- Bleaching; COD, Steady-state modeling: 
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1, INTRODUCTION 

Chemical pulps are bleached by the addition of 
chemicals chlorine dioxide or 
hypochlorite. The bleaching process removes lignin 
from the pulp, which reduces the yield of a pulp 
produced from a given initial quantity of wood. 
Bleached pulps generally give stable highk-brightness 
pulps. The concern about the environmental effects of 
chlorinated organics in bleach plant effluents is steadily 
increasing. In the interest of reducing environmental 
impact of pulp mill effluents, many companies are 
examining ways to modify effluent composition and 
reduce effluent volume. In the extreme, reduced effluent 
volume means a closed effluent-free pulp mill operation. 
With the chlorine-based bleaching processes currently 
in use, a closed mill is difficult to attain due to the 
corrosive behavior of the chlorine compounds when 
they are processed in recovery systems. Hence 
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operational analysis and optimization of process 
variables would enable in reduction of hazardous 
effluents. The available methods for analyzing the 
individual components responsible for pollution in other 
industries are also used for bleaching effluents. The most 
widely used characteristics include adsorbable organic 
halogen (AOX), total organic chloride (TOCL), 
biochemical oxygen demand (BOD), chemical oxygen 
demand (COD), total organic carbon (TOC), total 
dissolved solids (TDS), and color (Junna and Ruonala, 
11). 


Besides the respective proportions of lignin and 
bleaching chemical. there are a number of variables 
which affect the bleaching reactions and subsequently 
the pollution load. Most of these effects are common to 
all bleaching stages but their importance or function may 
vary from one stage to another. Mixing In any bleaching 
or brightening reaction, the homogeneity of the mixture 
composed of pulp(made up of liquid and solids) and 
bleaching agent (gas and or liquid). This involves proper 
pulp consistency contro! and an appropriate choice of the 
mixing technology. 


Pulp pre-treatment The behavior of a bleaching 
may be affected to varying degrees by the presence of 
pulping by products or metallic ions. Pulp washers 
solve in good part the potential problems associated 
to these elements. 


Temperature [t is generally agreed that temperature 
affects positively the quantity of lignin agent. 
Temperature generally enters the models as a 
nonlinear coefficient of a rate equation (Axegard et 
al. (3). 


pH The alkslinity or acidity of the solution as 
measured by its pH plays an important role in every 
bleaching stages (Daneault et al.(7); Axegard and 
Tormund (4); Gellerstedt et al. (8), Myers and 
Edwards (13)). It may act as a lever for modifying the 
reaction rate or act as an observer of some internal 
mechanisms of the reaction. The optimal pH varies 
from one stage to another. 


Pulp strength Preserving the strength properties of 
the bleached pulp is an issue for every bleaching 
sequence. Attacking pulp with a bleaching chemical 
may actually weaken its fibers if the chemical is inclined 
to react with the carbohydrates and break the cellulose 
chain. Selectivity, a term used to describe how the 
bleaching chemical prefers reacting with lignin instead 
of carbohydrates varies widely from one chemical to 
another and will thus be individually discussed. 


In the present study, a steady-state model of a four 
stage chlorination-extraction-hypochlorination 
Ihypochlorination 2 (CEHH) bleach plant has been 
developed for the calculating the pollution load in terms 
of COD. The model is verified and case studies for 
optimizing process variables are performed with the aim 
of decreasing the total amount of COD coming out from 
the four washers as shown in the Fig |. 


Il. MODELING OF THE BLEACHING STAGES 


The control of any process is generally preceded 
by a careful modeling of its behavior, When available, 
knowledge of the underlying principles (be they 
physical, chemical or mechanical) guide the investigator 
in choosing appropriate model structures for a given 
process. The stoichiometry and kinetics of many 
bleaching sub-processes have been investigated. 
Unfortunately, the mechanisms involved are extremely 
complex and rarely lend themselves to direct 
observation such that it is difficult to verify even well 
founded hypotheses. The degree of development of 
"bleaching models" is not uniform over the sequence of 
ing bleaching and brightening technologies, In the 
following study, the bleaching technologies have been 
loosely divided in "conventional" technologies chlorine, 
alkaline extraction, hypochlorite (CEHH), 
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Fig 1. Process flow diagram of a CEHH bleach plant 


The flow sheet of CEHH bleaching plant is shown in Fig. 
1. Each stage is composed of three unit operations as 
mixer, retention tower and washer. The process units and 
variables used to model each stage are shown in Fig. 2. 
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Fig 2: Process. units and variables of the bleach plant 


A. Mathematical Models of Unit Operations 


In a bleaching stage, mixers are used to mix the fibre 
suspensions with bleaching chemicals, steam or recycle 
streams. The mass balance equations in the mixer for the 
liquor, fibres, kappa number, and residual chemicals are 
presented in Table | based on the assumption of perfect 
mixing and quasi-steady-state, 


In the retention tower, the reaction between the bleaching 
chemicals and the fibres takes place. Wang (19), 
examined different models in order to represent the 
retention tower. Among the models that are tested, it was 
found that the sequence of CSTRS-PFR-CSTRS gives 
the best agreement with the experimental retention time 
distributions. Several researchers have reported two 
stage rate mechanisms for various stages of bleaching 
and if the concentration of reactants is very high in the 
beginning and fairly low in the end than it may be 
assumed that all retention towers are acting as single 
PFR. During further course of study all four towers are 
considered as PFR. The model equations are given in 
Table 1. 


The modelling of the drum washer is done by dividing it 
into three parts; as the dilution vat, drum filter and seal 
tank. In the dilution vat the liquor from the seal tank is 
mixed with the main pulp stream in order to dilute the 
pulp to the desired consistency. Therefore, it can be 
modelled as a steady-state mixer, The drum filter is used 
to remove the dissolved solids and unreacted chemicals 
from the pulp with wash water. In order to model the 
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drum filter, different models can be used. In this study, 
Norden efficiency approach (Norden, 14; Norden and 
Pekkanen, 15) is used in order to model the drum filter. 
In the model, the Norden efficiency factor is taken as 3 
for all the washers, In the seal tank the part of the liquor 
coming from the drum filter is recycled back into the 
dilution vat and the rest of it is sent to the sewer, which 
cause different environmental problems. Thus, the seal 
tank can be modelled as a splitter. The model equations 
for these unit operations are also given in Table 1 


B. Kinetic Models 


The kinetic models discussed in the literature try to 
describe and/or explain the chemical and physical 
interactions between the bleaching agent and the pulp 
fibers. In the modelling of the retention tower the most 
important parameter is the bleaching reaction kinetics. 
In the literature different kinetic models are proposed for 
the bleaching steps Ismail Dogan, et al (10). In the 
present study the process kinetic models are employed, 
using Kappa number for the delignification stages (C, 
E). However, in the brightening stages (H,, H,) light 
absorption coefficient is used, 


Most of the chlorination models are based on the kinetic 
study conducted by Ackert (1). In these equations the 
effect of consistency introduced through a rate constant 
multiplier, f, is equal to 1.0 for individual fibers (low 
consistency) and 0.52 for fiber mats (high consistencies 
where fibers contact each other).T is the temperature in 
s 0.5 (total initial lignin) and the 
is 0.3(total initial lignin), with non removable 
lignin = 0.2 (total initial lignin). These equations 
kinetically describe chlorine bleaching. The temperature 
dependence enters through the Arrhenius term for rate 
constant. The constant f;, is used to adjust the model for 
the low consistencies used for a conventional bleaching 
process as opposed to high consistency displacement 
bleaching. It should be noted that the pH dependence is 
not included in this model (Sklarewitz etal .18). 


Rage = KCl HL fas] 
Row =KLCh] Lond 


where 


fe 
K, = 0535 fexp, =| 
4 F 


(250 ) 
K ,= 0.0107 fexp! = | 
LT J 


Singh, S.V et al (17) conducted experiment on 
hardwood and proposed that the course of chlorination 
was governed by two stage rate equations: 


R =k (AKG, ol?) 


where c denote the concentration of available chlorine, c, 
the concentration of available chlorine during rapid 
phase of reaction and c, the concentration of available 
chlorine during slow phase of reaction . Singh, S.V et al 
(17) observed that the first term on right hand side of the 
equation (*) represented of the kinetics of the initial rapid 
reaction period up to 1 minute and the second term 
denoted the slow process and chlorination temperature 
did not effect the efficiency with which chlorine was 
consumed .this means that the slight variation in 
temperature during chlorination stage would not effect 
the pulp quality and its behavior in subsequent chemical 
stages of treatment. Chapnerkar (6) concluded the 
efficiency was also independent of both the initial 
chlorine concentration and consistency. Therefore, the 
kappa number after the chlorination as a function of the 
unbleached kappa number and the degree of chlorination 
a linear relationship is found in this case, which can 
represented by the equation (Axegard , P 2) 


Ke = K(unbleached) - 82 * consumed C/2 


* consumed C/, = % of 0.d. pulp / kappa number, Ke = 
kappa number after chlorination 


The kappa number after extraction stage could be 
calculated by using the rate equations given in table 2 (E- 
stage). Kappa number after the extraction stage (E) 
should be related with the inlet light absorption 
coefficient of the hypochlorination stage (H). Wang (19) 
found a correlation with the light absorption coefficient 
and the Kappa number of a bleached pulp by fitting the 
published experimental data. In our study, this 
correlation is used in order to link the Kappa number 
after the E stage and the light absorption coefficient 
entering the H stage as 


K,,, = 4.69 K, + .01 


After the extraction stage, the light absorption coefficient 
should be converted to brightness, because mills are 
using brightness for the quality of the bleached pulp, In 
order to convert the light absorption coefficient to 
brightness, the Kubelka-Munk equation was used 


K/S = (1-R,) /2R, 


where R, denotes the brightness, and S is the light 
‘scattering coefficient. The value of the scattering 
coefficient is used as 50 m/kg (Wang et al, 20). K is the 
light absorption coefficient of the pulp. 


TABLE | 


The Mathematical Models used in a Bleach Plant 


Liquors Fibers Kappa Number | Residual Chemical/Dissolved solids 
Mixer F,=F.+W, F(C,/l-C,)=FiG/I-C,) _[Ki=K, F,L,,=F,.L,+W4+M, 
PFR F,=K m= Oy dK/dt =-r dL /dt=-(Cy/1-C, wir 
Dilution Vat |F,=F,+W, Lows = (EL WME, 
Drum Filter |.............. Ls) = {Most RW)” (RW-1) 
(F, = 1)(Ligg,-Mog;)} When RW#1 
={Mpcx#(I/E,) )(Lins,Mosa) 
wine RW=1 
WP Wet cE, | ssisnesvivieuaracavens Mose = (WoMps2+F Lis .-F Lis) Wa 
SealTank [W,=W,+W, Mosa = Mos, 
TABLE 2 
‘The kinetic expressions used in the model 
Stage Reference Rate Expression Y 
Cc S.V.Singh et al (1975) R= -4#e =K, 6 +K,¢, 
where K,=Aexp # =0.85(min') at30° 
K,=Bexp (r]-226 10" (min %ofCL,)' at 30°C 
A=constant, B= constant, E = energy of activation , 
T= temperature in°Kk i 
z Hxegerd (197), Rue=KUOH } [say ) Rw =KLOH VU 
wherek,=1,96x10° ea 
K,=0.0103 exp =| 
Fun 4.75% 1o'exp{ 2358 ony, 
Chlorine Charge factor 
Axegard and Tormun R,~K[OCI]"L 0.688 
(1985) where K=7.9657 x 10"exp i“ 
TABLE 3 
‘The input data for simulation of CEHH bleach plant 
Cc E TH TH 
(1) Mill data Chlorination 
Flow rate (Kg/min) 4.4 - - > 
Exiration 
Flow rate (Kg/min) - 2.2 .25 25 
Concentration (%) - 8 8 8 
Hypo(1) as Cl, basis 
Flow rate (Kg/min) - - 2.4 - 
Concentration (%) - - 1 - 
. Hypo(2) as Cl, basis 
Flow rate (Kg/min) - - - 2.4 
Concentration (%) - - = 7 
Steam (Kg/min) - 12 - - 
Tower 
Temperature (OC) 25 65 45 45 
Residence Time , t (min) 70. 90 190 190 
Washer 
Wash liquor (kg/min) 1700 500 520 530 
Outlet consistency (%) 1 11.3 12 11.9 
(2) Unbleached pulp 
Pulp flowrate =111.14Kg/min| Consistency =3.0% 
Kappa number =25 Dissolved solids=0,001 % 


33 


34 


Model Equations for Pollution Load (COD) 


The equation for COD taken from Anjana (9) can be 
used for the mill data. The mill, already having ETP 
(effluent treatment plant) should consider issue of the 
impact of effluent load over and above the minimum 
carry over solids and difference of kappa number from 
floor kappa value of brown stock, Assuming the value of 
minimum carry over solids as mentioned in Anjana (9) 
to be a constant factor for all DF and N values, the Eq. 
for additional loads can be written as: 


COD = 2.4 * A Kappa Number + 0.61 (carry over solids) 


Where A kappa = K-K,, K,, being floor kappa number. 
K, depends on the type of fibers (HW, SW etc.) for a 
fully washed pulp, pulping condition, control 
measures, lignin content and type of lignin. Exact 
value of this parameter is generally not known. 


Il. MODEL VERIFICATION 


The present model is tested by comparing the simulation 
results with the literature mill data available in the 
Anjana (9), The input variables are taken from the 
literature and are tabulated in Table 3. The steady-state 
Kappa number, the residual chemical and the dissolved 
solids content of the filtrate from the washer for each 
stage are calculated and together with the literature data 
are presented in Table 4. A good agreement between the 
simulation and the mill data is achieved. 


IV. CASE STUDY 


The case study performed covers the optimization of the 
process variables. For each of the towers, the 
temperature, the concentration, the residence time and 
the consistency are the parameters investigated and their 
ranges are shown in Table 5, With in each range the 
variables were changed at constant intervals and the mid 
value corresponds to the steady state data of the bleach 
plant. Only for the chlorine flow rate the minimum 
allowable value is used. 


The results of the case study for the four stages are 
also given in Table 5, In this table, the first row indicates 
the steady-state literature data of the total amount of 
COD that was emerged from the four washers and the 
brightness at the end of the bleaching sequence. For the 
pulp to have an acceptable quality the brightness should 
be greater than 84% ISO. Effect of the flow rate of 
chlorine in the chlorination stage can not be studied due 
to lake of data. But it is observed that, in the mill, excess 
amount of chlorine is used in order to remove as much 
lignin as possible in the first stage. Therefore, the 
chlorine flow rate can be decreased without decreasing 
the brightness at the end of the sequence. Most of the 
coloring materials are removed in the first two stages The 
last two stages of the bleaching plant are mainly used to 
give the final brightness and to strength the pulp.. 
Therefore, the dissolved solids content of the filtrate is 
low in quantity for the last two stages. 


The results show that in each case the lowest 
allowable value of the parameters satisfying the 
brightness requirement should be used. Although for 
each stage the change of each parameter has a 
considerable effect, it has been compensated in the 
following stages. The pulp and paper industry consumes 
large volumes of water. An important measure in 
reducing the pollution charge is the reduction in water 
consumption and COD. 


In order to see the effect of each parameter, the 
optimized values of the parameters are simulated. The 
brightness at the end of the fourth stage is decreased to 
83.2%, still it is in the allowable range for the pulp 
quality. The bleaching processes mostly use molecular 
chlorine as the oxidizing agent to remove residual lignin 
from the cellulose fibres. The model that is proposed by 
Singh, S.V et al (17) is used in the simulation. The results 
obtained showed a 15% decrease in the total amount of 
dissolved solids witha brightness value of 85.37% ISO. 


TABLE 4 
Comparison of the mill data with the simulation data 
Stage Kappa no. or brightness Residual Chemical (%) COD (%) 
Simu. Litera. Simu. Report Simu Litre 
Cc 9.2 8-11 0.072 | 0.07 0.0528 0.0526 
E 3.03 2-7 0.261 | : 0.1763 0.1518 
HI 84.2 75-86 0.249 - 0.1141 0.1170 
H2 85.2 75-86 0,249 |. 0.1141 0.1170 
| 


TABLE 5 
Results of process parameter optimization in each stage 


Optimization CoD Brightness 
Parameter Range Value Red. (%) 
Steady state data - 85.37 
Chlorination 
Concentration (%) 6-8 1.465 85.37 
Temperature (°C) 20-40 20 0.089 85.37 
Consistency (%) 3-4 3 0.057 85.37 
Resid. time(min) 60-90, 60 0.046 | 85.37 
Extraction 
Concentration (%) 7-9 1 0.209 85.37 
Temperature (°C) 50-70 50 0.598 85.37 
Consistency (%) 10-11.5 10 0.510 85.37 
Resid. time(min) 90-110 90 0.423 85.37 
Hypochlorination(1) 
Concentration (%) 5-7 5 0.296 84.65 
Temperature (°C) 30-40 30 1.197 84.78 
Consistency (%) 3-5 3 0.243 84.69 
Resid. time(min) 170-210 170 0.176 84.81 
Hypochlorination(2) 
Concentration (%) 5-7 5 0.372 84.94 
Temperature (°C) 30-40 30 1.141 85.23 
Consistency (%) 3-5 3! 0.210 85.34 
Resid. time(min) 170-210 170 0.148 85.34 
Y. CONCLUSION E,, Norden efficiency factor (-) 

The steady-state model of a CEHH bleach plant was f, __rate constant multiplier 

asanenbled by linking ite a epseeon moses! for f. chlorine charge factor, % available C/, on 0.d. 

mixing, reaction and washing, Chemical oxygen 

demand (COD) is a widely used method for evaluating pulp/unbleached Kappa number 

the bleaching result, but as a collective measurement HE  contentofchiombphores exirieaeeRagoa 

variable it does not describe the actual compounds that bagi . 

cause the "loss" of bleaching chemicals. Studies have number in delignification stages (C and E) ml 

shown that many compounds contribute to COD load 0.1NKMn0,,, fibre or lightabsorption 

but ultimately most of them have no real effect on the . 

bleaching result. The applicability of the CEHH bleach coefficient in brightening stage 

lant model was checked by the literature data. A good § 
posueee between the edeet and the ifterature data (Hy, H) (m'/kg) 
was obtained. k rate constant of bleaching reaction 
Different case studies were performed in prder to F flowrate of liquor in pulp stream (kg/min) 

decrease the amount of COD coming out from the R, brightness (% ISO) 

washers. It was concluded that the effective process 

parameter for chlorination stage was concentration of | R — bleachingreactionrate, i.e. rate of Kappa 


bleach liquor where the temperature is more effective in 
Hypo stages. Extraction stage was the least effected by 
the above parameters. 


Nomenclature 


C, fibre consistency (% on slurry) 


number decrease ina delignification or rate 
of light absorption coefficient decrease ina 


brightening reaction 


RW. wash liquor ratio (-) 


S__ lightscattering coefficient (m‘/kg) 

T temperature (°K) 

t mean residence time in PFR (min) 

W — flow rate of liquor in non-pulp stream (kg/min) 

L dissolved solids or bleaching chemical content of 
the pulp stream (%) 

M dissolved solids or bleaching chemical content of 
the non-pulp stream (%) 

Subseripts 

C chlorination stage 

d drum filter stream 

DS dissolved solids 

E extraction stage 

H__ hypochlorination stage 

f fast 

s slow 

i into process unit 


J bleaching liquorcomponents C2, H, OH- and DS 


0 offprocess unit 
r recycle stream 
v_ dilution vat stream 
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